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Abstract. Given a network and a partition in n communities, we address the 
issues "how communities influence each other" and "when two given communities 
do communicate" . We prove that, for a small-world network, among communities, 
a simple superposition principle applies and each community plays the role of a 
microscopic spin governed by a sort of effective TAP (Thouless, Anderson and 
Palmer) equations. The relative susceptibilities derived from these equations 
calculated at finite or zero temperature (where the method provides an effective 
percolation theory) give us the answers to the above issues. As for the already 
studied case n = 1, these equations are exact in the paramagnetic regions 
(at T = this means below the percolation threshold) and provide effective 
approximations in the other regions. However, unlike the case n = 1, asymmetries 
among the communities may lead, via the TAP-like structure of the equations, to 
many metastable states whose number, in the case of negative short-cuts among 
the communities, may grow exponentially fast with n and glassy scenarios with 
a remarkable presence of abrupt jumps take place. Furthermore, as a byproduct, 
from the relative susceptibilities a natural and efficient method to detect the 
community structure of a generic network emerges. 
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1. Introduction 

Recently, in the network science, the issue to find the "optimal" community structure 
that should be present in a given random graph (a network) (£, F). C and T being the 
set of the vertices and bonds, respectively, has received much attention. The general 
idea behind the concept of community structure comes from the observation that in 
many situations real data show an intrinsic partition of the vertices of the graph into n 
groups, called communities, (£ = U"^j^£'^''' , F = U^^j.^j^F'^'^'^)), such that between any 
two communities there is a number of bonds that is relatively small if compared with 
the number of bonds present in each community. The partition(s) can be used to build 
a higher-level meta-network where the meta-nodes are now the communities (cells, 
proteins, groups of people, . . .) and play important roles in unveiling the functional 
organization inside the network. In order to detect the community structure of a 
given network, many methods have been proposed and special progresses have been 
made by mapping the problem for identifying community structures to optimization 
problems [2 [2], by looking for A;— clique sub-graphs 0], or by looking for clustering 
desynchronization [J. In general there is not a unique criterion. However, given 
a structure in communities, whatsoever be the method used, and assuming that 
the found partition (UjLj^£*^'\ U"<j.^]^r*^'''"')) represents sufficiently well the intrinsic 
community structure of the given network [5] , there is still left the fundamental issue 
about the true relationships among these communities. Under which conditions, and 
how much two given communities, also in the presence of other communities, exchange 
information, how they influence each other, positively or negatively, what is the typical 
state of a single community, what is the expected behavior for n large, etc... are all 
issues that cannot be addressed by simply using the above methods to detect the 
community structure. In fact, all these methods, with the exception of Refs. [1] and 
[3], are essentially based only on a topological (and, in most cases, local) analysis of the 
network. To uncover the real communication among the communities we have to pose 
over the graph (£, F) a minimal model in which the vertices assume at least two states, 
and analyze their correlations. Confining the problem to the equilibrium case we have 
hence to use the Gibbs-Boltzmann statistical mechanics. Given a community structure 
of size N, if each node is associated with a spin cr^, i = 1, . . . , iV, then each meta-node 
is associated with the meta-spin s*-'' = Eie£(0 ^ ~ 1, ■ ■ ■ , where the 

sums run only over the nodes of the ^-th community. We show then that the average 
magnetizations m*^'' of the meta-spins obey special effective equations as if they were 
microscopic spins immersed in a ferro or glassy material. From these equations it 
will be then simple to derive the relative susceptibilities x*-'"'^-', I, k = 1, . . . ,n among 
the communities revealing in an unambiguous way the communication and correlation 
properties of the meta-network. 

2. A minimal model 

If a Gibbs-Boltzmann distribution, exp(— /3iJ), has been assumed, H being some 
effective Hamiltonian, one can obtain the (adimensional) couplings pjf'^'' from the 
data of the given graph by isolating the two vertices i,j from all the others, and by 
measuring the correlation function of the obtained isolated dimer {uiOj)' , where (■)' 
stands for the Gibbs-Boltzmann average of the isolated dimer. The general problem 
is actually more complicated due to the presence of two sources of disorder since 
both the set of bonds F, and the single couplings {jff^}, may change with time. 
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Assuming that the time scale over which this changes take place is much larger than 
that of the thermal vibrations of the spins, we have then to consider a disordered Ising 
model with quenched disorder. Here we specialize this general problem to the case of 
Poissonian disorder of the graph, while we leave the disorder of the couplings arbitrary. 
We formulate the problem in terms of Ising models on generic small-world graphs 
[B]: given an arbitrary graph {£,Oi^o)i the pure graph, and an associated community 

structure {^Y=i^o'^ ^^?<k=i^o''^^)^ which each community has an arbitrary size, we 
consider a generic Ising Hamiltonian Hq defined on this non random structure, the 
pure model, characterized by arbitrary couplings Jg'"'"'' , and, between any two sites 
i,j, with i £ jC^'''> and j S jC^''\ we add some random connections (short-cuts) having 
average connectivities c'-'''^^ (more precisely, if |i3g'^'' | is the size of the k-th community, 

(I k) (I k) 

we introduce the directed random variables q j taking the values ] =0,1 with 

probabilities 1 — c'^''''^'^ /IC^'^ \ and c^'^'^-'/l'^o'^'' I? respectively). Then, along with these 
random connections, random couplings jC^'^) with a quenched disorder arc imposed, 
and we study the corresponding random Ising model, the random model, having 
Hamiltonian H (see Fig. 1). In a compact way H is therefore given by 

n 

H {{a.,}ti) = Ho {{a.}l,) + AH ({ajf^,} - '^''^ E (1) 

In Eq. (P) Ho is the non random part of the Hamiltonian having non random couplings 
Jq''*^'' (typically short-range couplings but not necessarily) whereas AH is the random 

(I k) (I k) fi\ 

part of the Hamiltonian involving only the long-range couplings c] j J^j , finally ' 




Figure 1. (Color on-line) An example with n = 2, l/^g"*"'! = 6, |£q^^| = 4, 

IFq ' '^I =5, IFq ' 'I =4, Irjj ''^1 = 1. Continuous lines represent short-range 
(I k) 

couplings Jg ' of the pure model, whereas dashed and dot-dashed lines represent 
the additional random couplings , with i € £g ' and j G ■Cg'°' , and are present 

only in the random model (panel c)). Panel a) and b) differ for the position of 
f 1 2) 

one single coupling Jg ' , nevertheless the example in a) , compared with the case 
in b), at least at finite temperature (where the lengths of paths are important), 
benefits clearly of a better communication (higher betweenness) . 
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is an arbitrary external field acting only on the Z-tli community. 
3. Equations of the meta-network 

In [7] we established a new general method to analyze critical phenomena on the 
small-world models represented by Eq. ([T]) for the case in which there is one single 
community, n = 1: given any arbitrary initial graph {Co, Fq), if we increase the average 
connectivity by c through the addition of Nc bonds randomly spread over the initial 
graph, under the only condition that be c > 0, we have found an effective field theory 
that generalizes the Curie- Weiss mean-field theory via the equation 

m(^) =mo(/3jf\/3j(^)m(^)) (2) 

and that is able to take into account both the infinite and finite dimensionality 
simultaneously present in small-world models. In Eq. ^ mQ(f3Jo] (3h) represents 
the magnetization of the pure model, i.e., without short-cuts, supposed known as a 
function of the short-range coupling Jq and an arbitrary external field h, whereas the 
symbol E stands for the ferro-like solution, E = _F, or the spin glass-like solution, 
E ~ SG, and Jq^' and J*-^-* are effective couplings. Here we generalize this result to 
the present case of n communities of arbitrary sizes and interactions; short-range and 
long-range (or short-cuts) couplings. We show that, among the communities, a natural 
superposition principle applies and we find that the n order parameters, F or SG like, 
obey a system of equations which, a part from the absence of the Onsager's reaction 
term [8], can be seen as an n x n effective system of TAP (Thouless, Anderson and 
Palmer) equations [HI in which each community plays the role of a single "microscopic" - 
spin through its own order parameter, nS^''^\ I = 1, . . . , n flOj . In particular, in the 
simpler case in which there are no short-range couplings among different communities 
[J^'^^ ~ Q ior I ^ k) these self-consistent equations take the form 

n 

fc=i 

where ttiq ^(/3Jq'^; (3h^^^) is the magnetization of the pure model for the l-ih community, 
and the effective couplings jt^-'-'^), jiSG;i,k}^ Jq^'^^ and Jq^^''-* are given by 

pjiF;i,k) dg ^(;,,) J ^^(;^,) tanh(/3 J) , 
^j(SG;/,fc) dg J ^^(,^,) tanh2(/3J), 

4\ «^') tanh-i(tanh2(^jW)). 

In the above definitions dfi^''''''' (J) stands for the probability distribution of the long- 
range coupling disorder between the ^-th and k-th community. It is easy to see that, 
if the size of the communities is parameterized as |>Cq'*| = a'-'^iV, with Ym=i '^^^^ ~ 1' 
then the connectivities must satisfy the balance equation a'-'^c'''''^'' — a'^'^^ cJ'^'^\ so 
that, when a*^'^ ^ a^''\ the effective couplings and j(SG;;,fc) ^^^^ j-^q^ symmetric 

even if the random couplings J^'j''^ are symmetric. Though being not complex, the 
rigorous derivation of Eqs. ^ is quite lengthy. We refer the interested reader to 
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|10j . However, if we make the natural assumption that the effective couphngs satisfy 
a Hnear superposition principle, we then see that Eqs. ^ are immediately derived 
from Eq. 

As for one single community, Eqs. ^ are exact in the paramagnetic region (P) 
while provide an effective approximation in the other regions. More precisely, off the 
P region, for unfrustrated disorders Eqs. ^ are exact up to 0(l/c^'''^^) terms and 
become exact also in the limit c'''*^' 0+, whereas for frustrated disorders Eqs. ^ 
in general give only a qualitative effective description of the order parameters. We 
stress however that in both the cases the critical surfaces derived from our effective 
Eqs. are exact (notice in contrast that the pure naive mean field equations give 
a wrong critical surface) whenever, for any ^ = 1, . . . , rt, it is c'^'''-' > or, if for some 
I, c^'''-* = 0, there is at least a chain of, say h, connected communities such that 
(.{i,h) > Q^c''hM > 0, . . . , c^''"''') > 0. We remind the reader that the remarkable 
progresses reached in the context of mean field theory for disordered models, as in 
the case of the ordinary TAP equations (or the Viana-Bray model [H]), concern only 
models where short loops in the graph are absent, while in our approach the graph 
{Co,Tq) may be an arbitrary lattice, regular or not, and having loops of any length. 
The possibility to improve the theory off the P region is a formidable task to be 
investigated in the future, however, as we show below, it should be already clear the 
very interesting phase-transition scenario emerging from Eq. ([3]). 

In [7] {n = 1) we established the general scenario of the critical behavior coming 
from Eq. stressing the differences between the cases Jo > and Jq < 0, the 

former being able to give only second-order phase transitions with classical critical 
exponents, whereas the latter being able to give rise, for a sufficiently large connectivity 
c, to multicritical points in principle also with first-order phase transitions. The same 
scenario essentially takes place also for n > 2 provided that the Jq's and the J's 
be almost the same for all the communities, otherwise many other situations are 
possible. In particular, unlike the case n = 1, relative antiferromagnetism between 
two communities I and k is possible as soon as the Jj''^-'^'' have negative averages, while 
internal antiferromagnetism inside a single community, say the l-th one, due to the 
presence of negative couplings J^^ < 0, is never possible as soon as disorder is present 
[12j . Less intuitive and quite interestingly, if we try to connect randomly with some 
added connectivity c'^'^'^' the l-th community having inside only positive couplings 
("good") to the fc-th community having inside only negative couplings ("bad"), not 
only the bad community gains a non zero order, but even the already good community 
gets an improved order. In Fig. 2 we report an example. However, with respect to 
the case n = 1, another peculiar feature to take into account is the presence of many 
metastable states. In fact, this is a general mechanism of the TAP-like structure of 
the equations: as we consider systems with an increasing number of communities, 
the number of metastable states grows with n and may grow exponentially fast in 
the case of negative short-cuts. A metastable state can be made virtually stable (or, 
more precisely leading) by forcing the system with appropriate initial conditions, by 
fast cooling, or by means of suitable external fields. As a result, with respect to 
variations of the several parameters of the model (couplings, connectivities, sizes of 
the communities) , the presence of metastable states may lead itself to first-order phase 
transitions even when the Jo's are all non negative. This general mechanism has been 
already studied in the simplest version of these models, namely the n = 2 Curie- Weiss 
model {Jq''''' = and c^'''"'^ — > oo), where a first-order phase transition was observed 



Small-world of communities: communication and correlation of the meta-network 6 

1 




(F;l), (1.2) 

m (c 



(1 2) 

L(c* ' =0) 

(F;l), (1.2 

m (c 



Figure 2. (Color online) Solutions of Eqs. ([Sj with n = 2 in the case Jq ' = (a 
generalized Viana-Bray model) with long-range couplings = 1, = —1 

and = 1 (j'^'^' = —1 also leads to the same plot), and connectivities 

c(i>i) = c'^'^' = 2, and c'-"^'^) = or c'-*^'^' = 1. L stands for the free energy term 
associated to each solution. 




Figure 3. (Color online) Solutions of Eqs. lO with n = 3 in the case 

and c^''*) — > oo (a generalized Curio-Weiss model), with jC'') = S/c^'''^ and 

j(i.fc) = — 0.5/c('''^\ for I ^ k. We plot also the associated free energy term L. 



to be tuned by the relative sizes of the two communities and by the external fields 
|13j : moreover, first-order phase transitions have been observed in simulations of a two 
dimensional small- world model with directed shortcuts [14j . In particular, in system 
of many communities, n 3> 1, a remarkable and natural presence of first-order phase 
transitions (tuned by the several parameters) is expected which, if the J's or the Jo's 
arc negative, reflects on the fact that the communities, at sufficiently low temperature, 
behave as spins in an effective glassy state [TS1[TB|. Fig. 3 concerns the Curie- Weiss 
case with n = 3; as is evident, even for small n, the number of mctastable states is 
rapidly growing. 
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4. Communication properties of the meta-network 

In general, how much two communities influence each other, is encoded in the matrix 
^(i,fc) '^S/ the (adimensional) susceptibility of the random model 
which tells us how the l-th community reacts to a small variation occurred only 
(initially) in the fc-th community. From Eq. ([3]), or from its most general form that 
includes also non zero short-range couplings among different communities, we have 

X(^) = (l-X0-/3J(^))"'-X0, (4) 
where we have introduced the matrix of the effective long-range couplings f3J^^'^'^\ 

~(l k) . . 

and Xq , the adimensional susceptibility of the pure model. Note that in the case 
4^''"^ = for ; ^ fc, Xo is a diagonal matrix whereas x^^^ is not. By looking for the 
points where x'^^^ becomes singular, from Eq. Q we see immediately that the critical 
surface Pc where the second-order transition takes place is solution of the following 
exact equation 

dct(l-xo-/3J^^^) =0. (5) 

By sending (3 oo, the theory can be projected in particular at zero temperature 
where, for positive couplings, a natural effective percolation theory arises. Then, by 
using lim^^+oc /3 J^^^ = c, in the limit (3 ^ oo Eq. ^ becomes 

£ = {l-£o cy' So, (6) 
where we have introduced 



{l,k) def 




^ii„^ir'({rf^'''}^o). (7) 

Eq. ([6|) tells us how £ changes as we vary c, being an exact equation as c belongs to 
the P region, i.e., below the percolation threshold Cc, solution of the following exact 
equation [J 

det {l-£o-c) = 0. (8) 

At T = there is no thermal dissipation and it is easier to analyze the 
communication properties. Given the arbitrary pure graph [Cq, Fq), and a community 
structure assignment which spHts the set of the bonds Fq in n{n — l)/2 sets, 

(Ik). 

To = Ui^feFg , in the pure graph the communities I and k communicate if and only if 
£q'''^ _^ whereas in the random graph if and only if fC^'^' ^ 0. We can understand 
the communication process by observing that the characteristic time to exchange 
a unit of information between the two communities / and k in the pure model grows 

as tg'^'' oc (^Sq''^^^ and, similarly, for the random model as i^'"'"'^ oc (f^'''^)) ^ . From 

Eq. ([6|) we see that in the pure model, if, for I ^ k, £^^'^^ ~ 0, the two communities 
/ and k cannot communicate (ip '^"^ oo), but for any arbitrary small c'''*^' > they 

X In Eq. Q it is understood that we are considering only graplis {Cq, Fq) for which the pure model 
has no critical temperature even for T = 0. In fact, if this is not the case, as occurs for instance if 
{Cq,Tq) is a d(j dimensional lattice with do > 1, the P region is shrunk to the single trivial point 
^(i,fc) ^ Q jj^j^j take effectively S^'^^ = +oo |10) . More precisely, if /3cO < cxd or even if f3c0 = 0, 

as occurs in scale free networks with exponent 7 < 3 1171 . one should use Eqs. (O keeping A'^ finite. 
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communicate and the characteristic time decays with c approximately as (it is easy 
to see that for any I is always Sq'''^ > 1) H 

(^^(M)^(W^(fc,/c)^^(^2)^-\ (9) 

whereas, at higher order in c, Eq. ^ takes into account that the communities / and 
k can communicate also indirectly via chains of other intermediate communities. In 
general, if the pure graph has an own dimension do (possibly fractal) sufficiently 
high, do > 1, one has fg''^^ ~> oo for N —^ oo, so that we have tg''"^ — > 
and then also tC^*^) — > 0; i.e., the communities communicate instantaneously (they 
percolate). However, in the random model, even if fg''*^-* is finite, when c approaches 
the percolation threshold surface Cc, then we have t^'^'^-' ^ 0, a peculiar feature which 
is possible only in the random model. We stress again that the graph (£g,rg) is 
completely arbitrary. So, for example, for n = 1, it is easy to check the consistency 
of our effective percolation theory (i.e., we recover the same critical surface) with the 
recent model introduced in |18| where the classical Erdos-Rcnyi random graph |19| 
is generalized to include a finite clustering. However, whereas in [T^ the percolation 
analysis is performed by ab initio calculations starting directly from graph theory 
elements (at least this seems in principle possible for graphs having regular loops), 
our effective percolation theory requires to perform a simulated annealing toward 
T = of the non random model defined over the pure graph (£g,ro) and immersed 
in a small external field. At each small but finite T we make the simulation, then, 
once the relevant quantities like xa or fno are obtained for small T, the percolation 
properties of the random graph (C, T) are easily calculated. Therefore, in our effective 
percolation theory, in Eq. ^ the matrix £o represents an input data. In general, it 
can be sampled efficiently by simple simulated annealing procedures by using Eq. ([7]) 
since the problem is mapped to an unfrustrated Ising model (/3Jq''*'^ > 0). 

The matrix Sq (or, more in general, at finite T the matrix x), leads also to a 
natural new criterion to detect community structures: given an hypothetical value 
n, by a suitable generalization of the modularity introduced by Girvan and Newman 
P] which makes use of £o rather than that of the adjacency matrix, we can define 
a "measure" which takes into account paths of arbitrary length rather than links, 
and the resulting community structure coincides with that partition such that the 
communication among all the communities is minimal [lOj . In the algorithm proposed 
in [2], given the network, one removes the link having the highest betweenness, where 
the concept of betweenness, a measure of the centrality of the given link, can be 
given in several ways. In particular the betweenness of a link can be defined as the 
number of geodesic paths (the shortest path connecting two sites) passing through 
it. After deleting many times the links having the highest betweenness, a partition of 
the original network in communities can be obtained and a measure of the quality of 
assignment in communities is given as 



(10) 



where e^'''^'' is the fraction of all bonds connecting the commimities I and fc, and aj'' is 

defined as a^^ '''= ^f^e-^i^\ The term (a^'^)^ in Eq. pU)) represents the expected 
fraction of bonds falling inside the community I when their ends are connected 

§ Analogous relations hold also at finite T, provided the source of the signal be sufficiently slow. 
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randomly. Thanks to the presence of the term (a^'')^ in Eq. (fTO|) . Qi gives measure 
when one considers the trivial case in which Tq is a single community {n = 1), 
and partitions that maximize Qi correspond to best community structures. We can 
consider however other similar measures that takes into account not only bonds, but 
also, for example, paths of two consecutive bonds. In general we can define 

, 2" 



E 



— . . 

where now e^y^^^ is the fraction of all paths of length not greater than h connecting 

the two communities I and fc, and cif^ 'Ylik^h^^ ■ Again we have that its square 
represents the expected fraction of paths of length not greater than h having both 
ends inside the community / when they are connected randomly, and makes the 
measures (fTTj) non trivial. When h oo, the matrix e is proportional to the matrix 
£o. It is important to note that, at T = 0, the algorithm we propose to detect 
a community structure coincides with that of Newman and Girvan for the case of 
geodesic betweenness (this can be seen from the combinatorial meaning of the matrices 
fg'''^'' or [TD]), but the measure associated to the found partition is given by using 
Qoo, and not Qi: essentially Qi measures the relative lack of links among different 
communities, while Qoo measures the relative lack of communication among different 
communities. 

We point out that, in the pure model, having some bonds between the l-th and 

..(Ik) 

k-th communities does not guarantee that the condition £q > be satisfied. In 
fact, it is not difficult to see that to have Sq''^'^ > it is necessary that the number 
of paths between the l-th. and the k-th communities be at least of order N [TU] . Note 
also that such a requirement does not exclude the possibility that even a single bond 
between the two communities be enough, provided that through this bond passes 
at least 0{N) paths (high betweenness, or centrality; see Fig. 1). It should be 
then clear that measures based on a local analysis, and on elementary use of the 
adjacency matrix, cannot capture the real communication properties. We illustrate 
now in a simple example, analytically feasible in our approach, how remarkable can 
be the differences in taking into account just links or, more properly, paths of any 
length. Such differences become actually much more important and interesting at 
finite temperature, where the length of each single path affects the susceptibility, while 
at T = the length of a path does not play any role. However, for simplicity in the 
following we will consider only the case T = 0. Suppose we have found n Erdos-Renyi 
sets, having intra averages connectivities c^'"'' and connected with each other randomly 
with inter averages connectivities c'-'''^^, I ^ k. In this case, we have Xo = 1, therefore, 
from Eqs. © and ([7]), we get immediately £ = (1 — c)~^, which is very different from 
the adjacency matrix c. This can be well understood by observing that, below the 
percolation threshold surface given by Eq. ([HI) , it holds £ = (1 — c)^^ = 1 + c + c^..., 
where each term in the sum takes into account the presence of paths of length 0, 1, 
2,. . . respectively. Let us apply this result to the following six examples with n = 2 
and n = 3 communities: 

0.35 0.3 \ / 0.4 0.3 \ / 0.6 0.3 \ / 0.7 0.3 \ 



I 0.3 0.35 j I 0.3 0.3 j \ 0.3 0.1 / \ 0.3 0.0 
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and 




0.4 0.1 0.1 \ / 0.6 0.1 0.1 \ 

0.1 0.2 0.1 C6 = 0.1 0.1 0.1 
0.1 0.1 0.1 / \ 0.1 0.1 0.0 / 



In all the above examples the matrix c is always normalized, so that we have always 
e = c. Concerning the matrix Coo, it coincides with the matrix £ = {1 — c)~^ to be 
normalized to 1. The results corresponding to the above six examples, together with 
the modularities Q and Qoo, and the maximum eigenvalue A of the matrix c, give, 
respectively: 

Q = 0.2, Qoo = 0.312, A = 0.650 
Q = 0.195, Qoo = 0.310, A = 0.654 
Q = 0.007, Qoo = 0.264, A = 0.740 
Q = -0.045, Qoo = 0.217, A = 0.810 
0.042 \ 

0.032 Q = 0.320, Qoo = 0.463, A = 0.487 
0.221 / 

0.049 \ 

0.025 Q = 0.160, Qoo = 0.418, A = 0.640. 
0.172 / 

From these examples we observe the following: i) as expected Q and Qoo tend to be 
"parallel" measures, but they are quite far from being proportional; ii) from the point 
of view of communication, these six networks do not present very strong differences (we 
expect stronger differences in other kinds of networks embedded in some geometry; 
as the case depicted in Fig. [T|); Hi) as we pass from the case 1 to the case 4, and 
similarly from the case 5 to the case 6, i.e., toward more and more asymmetric cases, 
both Q and Qoo take lower and lower values, meaning that the partition has a poor 
meaning (Q small), or that the global communication among the communities is higher 
(Qoo small). This latter observation is less trivial and deserves attention: as a rule, 
asymmetric situations benefit of a better communication. Consider in particular the 
cases 5 and 6 and compare the matrix c (or e) with the matrix Coo: despite e^^'^^ and 
g(i.3) equal, we have that e^''^'' is almost twice 6^^^- This is due to the fact that, 
between the second and third community, the latter benefits of a stronger asymmetry 
with respect to the first community, which in turn has the largest density of intra 
bonds. For the same reason, in the case 6, we note also that the matrix element e^'"^'' 
is of the same order of magnitude of e^'^^ , despite being e^^''^^ = 0. 



_ / 0.406 0.094 \ 

~ 0.094 0.406 ) 

_ ( 0.437 0.094 \ 

~ 0.094 0.375 ) 

_ ( 0.562 0.094 \ 

eoo;3 - Q Q j 

_ f 0.625 0.094 \ 

~ 0.094 0.187 ) 

and 

/ 0.333 0.023 

eoo;5 = 0.023 0.249 

\ 0.042 0.032 

/ 0.436 0.027 

eoo;6 = 0.027 0.191 

\ 0.049 0.025 
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5. Conclusions 

Real- world networks present an intrinsic partition in communities. However, despite 
important progresses, the absence of universality of the many proposed techniques, 
makes the community detection an "Art" rather than a solid science [5D]. One weak 
point of these techniques lies on the fact that in most cases only the topology (i.e., 
the structure of the graph) - and often only a local topology - is taken into account 
and correlations never introduced. In particular, the real communication properties 
among the communities cannot rely on a local analysis. On the other hand, starting 
from real-data it is possible to define in an unambiguous way a minimal model, a 
disordered Ising model, able to take into account all the correlations, short- and long- 
range like, present on the given network. We then discover that, whatever be the 
given community structure, the exact relationship of the meta-nodes are regulated by 
a quite universal form of effective TAP equations which, through F- and SG-like order 
parameters and then the matrix x of the relative susceptibilities (in general completely 
different from the adjacency matrix), give rise to a reach variety of configuration 
and communication scenarios which are analytically and/or numerically feasible. In 
particular, by simulated annealing procedures applied to the non random model 
(c = 0), Eqs. (O-Q allow us to analyze how the percolation and communication 
properties of the system changes when it becomes random (c 7^ 0). We can also 
use the same technique to analyze the non random graphs itself for both studying 
the communication properties or its community structure. We stress that there is no 
limitation in the choice of the non random graph [Cq, Tq) and in particular, as occurs 
in real-world networks, it can contain loops of any length where therefore traditional 
mean-field theories developed for disordered systems (that suppose a tree-like structure 
of the graph) cannot be applied. Our approach is exact in the P region (where however, 
we recall, correlations of two spins neighbors of a same spin are not zero due to the 
presence of short loops) and on its boundaries (the critical surface). The possibility 
to improve the method off the P region is under investigation. We finally anticipate 
here that it is possible to generalize all the results to the cases of small-world models 
with free-scale intra- and inter-connections among communities [21] . 
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